gmt grdimage dem.grd -JM24 -B0.25i0.25i -Xc -Yc -Cmby.cpt -V > colored.ps
gmt grdgradient dem.grd -Ghillshade.grd -A345 -Ne0.6 -V
gmt grdhisteq hillshade.grd -Ghillshade-hist.grd -N -V
gmt grdinfo hillshade-hist.grd
gmt grdmath hillshade-hist.grd 5 DIV = hillshade-int.grd
gmt grdimage dem.grd -JM24 -B0.25i0.25i -Xc -Yc -Ihillshade-int.grd -Cmby.cpt -V > shaded.ps
S1A_IW_SLC__1SDV_20190806T224456_20190806T224523_028455_033731_BD53 S1A_IW_SLC__1SDV_20200309T224455_20200309T224522_031605_03A45F_C35A
# Executed
time p2p_S1_TOPS_Frame.csh \
S1A_IW_SLC__1SDV_20190806T224456_20190806T224523_028455_033731_BD53.SAFE \
S1A_OPER_AUX_POEORB_OPOD_20190826T120644_V20190805T225942_20190807T005942.EOF \
S1A_IW_SLC__1SDV_20200309T224455_20200309T224522_031605_03A45F_C35A.SAFE \
S1A_OPER_AUX_POEORB_OPOD_20200329T120846_V20200308T225942_20200310T005942.EOF \
config.s1a.txt vv 1
# real 472m6.076s
# user 528m19.833s
# sys 1m21.343s
library(rgdal)
library(raster)
setwd(paste0(
'/media/jose/BUDOCS/gmtsar_course_images/processed/',
'Santo_Domingo_all_subswaths_processed_20190806T224456_20200309T224455',
'_based_on_topex_gmtsar_example_S1A_SLC_TOPS_LA/merge/')
f <- raster(readGDAL('phasefilt_ll.grd'))
f <- raster(readGDAL('corr_ll.grd'))
c[c<0.2] <- NA
fm <- f*(c>=0.2)
extent(fm)[1:2] <- extent(fm)[1:2] - 360
writeRaster(fm, 'phasefilt_mask_ll.tif')
S1A_IW_SLC__1SDV_20160318T224427_20160318T224457_010430_00F779_D477 S1A_IW_SLC__1SDV_20180119T224443_20180119T224510_020230_022857_A730
# Executed
mkdir sd1618
cd sd1618
mkdir topo
mkdir raw
cd raw/
ln -s ../orig/*.tiff .
ln -s ../orig/*.xml .
ln -s ../orig/*.EOF .
ln -s ../topo/dem.grd .
time align_tops.csh \
s1a-iw1-slc-vv-20160318t224427-20160318t224455-010430-00f779-004 \
S1A_OPER_AUX_POEORB_OPOD_20160407T121449_V20160317T225943_20160319T005943.EOF \
s1a-iw1-slc-vv-20180119t224445-20180119t224510-020230-022857-004 \
S1A_OPER_AUX_POEORB_OPOD_20180208T120528_V20180118T225942_20180120T005942.EOF \
dem.grd 2>&1 | bash -c 'trap "" INT; awk "{print strftime(\"%c - \") \$0}"' | tee log
# real 2m15.453s
# user 2m43.789s
# sys 0m2.241s
cd ..
rm -r F1/raw
mkdir F1
cd F1
ln -s ../config.s1a.txt .
mkdir raw
cd raw
ln -s ../../raw/*F1* .
cd ..
mkdir topo
cd topo
ln -s ../../topo/dem.grd .
cd ../..
cd F1
time p2p_S1_TOPS.csh S1_20160318_224427_F1 S1_20180119_224445_F1 config.s1a.txt 2>&1 \
| bash -c 'trap "" INT; awk "{print strftime(\"%c - \") \$0}"' | tee log
# real 26m2.930s
# user 26m12.301s
# sys 0m7.845s
LOS map generated using the config file from Sentinel-1 TOPS Los Angeles example, which includes interpolation in SNAPHU. I avoided using the Greece example config file. See also this and this for more on horizontal-vertical decomposition
| Fecha | Parte | Precio |
|---|---|---|
| 7/jul/2019 | Raspberry Pi 3 B+ | US$38 |
| 8/jul/2019 | u-blox NEO-M8T | US$75 |
| 16/jul/2019 | Pantalla táctil Waveshare 480x320 | US$33 |
| 17/ago/2019 | Atenna TOPGNSS (sustituyó antena u-blox) | US$53 |
| 18/ago/2019 | Impresión 3D de caja, repo GitHub de Taroz | US$37 |
| 17/ago/2019 | 2 Baterías 8800mAh NP-F970 (Sony) con cargador | US$38 |
| 21/ago/2019 | Adaptador para batería, conversor a ~5V, cables | US$28 |
| 13/sep/2019 | Trípode | US$45 |
| 25/oct/2019 | Palo 2 m con nivel de burbuja | US$65 |
| Total sin ensamblar | US$161 | |
| Total ensamblada | US$302 | |
| Total ensamblada con trípode | US$347 | |
| Total ensamblada con trípode y palo | US$412 |